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ABSTRACT 

We study the angular clustering of ~ 6 x 10® NVSS sources on scales ^ 50h~^ Mpc in the context 
of the ACDM scenario. The analysis partially relies on the redshift distribution of 131 radio galaxies, 
inferred from the Hercules and CENSORS survey, and an empirical fit to the stellar to halo mass 
(SHM) relation. For redshifts z iS 0.7, the fraction of radio activity versus stellar mass evolves as 
where ao = 2.529±0.184 and ai = 1.854lg/g®. The estimate on is largely driven 


by the results of Best et al. (2005), while the constraint on ai is new. We derive a biasing factor 
h{z — 0.5) = 2.093('(g;igg between radio galaxies and the underlying mass. The function b{z) = 
0.33z^ -|-0.85z-|-1.6 fits well the redshift dependence. We also provide convenient parametric forms for 
the redshift dependent radio luminosity function, which are consistent with the redshift distribution 
and the NVSS source count versus flux. 

Subject headings: Cosmology: large scale structure of the Universe, dark matter 


1. INTRODUCTION 

Radio activity is a marker for energetic processes re¬ 
lated to star formation and accretion of gas on super 
massive black holes (SMBHs) at the nuclei of galaxies. 
Active Galactic Nuclei (AGN) are associated with cycles 
of powerful jets generated by accretion disks on SMBH. 
These jets make a sizable contribution to the energy bud¬ 
get at moderate redshifts and can play an important role 


et al. 

2006 Croton et al.||2006p and e 

ven ottsetting gas 
s of galaxies (e.g. 

cooling in massive groups and cluster 

Boehringer et al.||1993[ |Carilli et al.| 
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is of great interest 
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ment, whether galaxies are isolated or in dense regions, 
may also have an important impact on the emergence of 
radio activity. Galaxy-galaxy interactions tend to boost 


accre t ion of gas on the cen tral SMBH (e.g. Best et al. 
2005| |Pasquali et al.||2009|), eventually triggering AGN 
activity. One approach is to directly as sess the observed 
galaxy env i ronment of radio AGNs (e.g. Prestage fc Pea- 


cock||198^ IHill fc Lilly||1991[ IMcLure fe Dunlo 


Best et al.||2005| |Worpel et al.||2013p, but this is typ- 

ically limited by small sample uncertainties. Alterna¬ 
tively, the clustering of radio galaxies quantified through 
two-point statistics (e.g. correlation functions), could 
constrain the way these galaxies are related to the un- 
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are currently unsuitable for a robust determination of 
two point statistics at moderate redshifts. Here we re- 
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visit the angular power spectrum computed from radio 


sources in the N RAO VLA Sky Survey (NVSS) (Con- 


don et ar][l998|). The survey provides angular positions 

nr 


and 1.4 GHz huxes, S, for about 80% of the sky and is 
currently the most suitable survey for angular angular 
clustering analysis. 

Two ingredients are needed for extracting information 
from angular clustering statistics within a cosmologi¬ 
cal scenario, e.g. the ACDM. The first is a model for 
the mean redshift distribution, N{z), of NVSS sources 
(galaxies). We base our modeling of iV(z) on two catalogs 
of radio sources with redshifts but small angular cover¬ 
age. T hese are the Combined EIS-NVSS S urvey (CEN¬ 
SORS) (Best^^^^^^Rigb^et^ 2011) and the Her¬ 
cules survey (jWaddington et al.|2UUl ) which are (jointly) 
complete down to a hux limit ot 7.2mJy. The second 
ingredient is the biasing relation between the spatial dis¬ 
tribution of radio galaxies and the underlying mass den¬ 
sity field. Our modeling of this relation adopts a para¬ 
metric form for the frequency of radio activity versus 
stellar mass and an observationally motivated stellar to 
halo mass (SHM) ratio from the literature. Well estab¬ 
lished halo biasing is then incorporated to yield a pre¬ 
scription for the biasing relation of radio galaxies. The 
angular power spectrum of radio galaxies as predicted 
in the ACDM is then contrasted with the clustering of 
NVSS galaxies in order to put constraints the biasing re¬ 
lation and the prevalence of radio activity as a function 
of the stellar mass and redshift. 

The outline of the paper is as follows. In 0 we re¬ 
view the definition of the angular power spectrum, C;, 
and how it is estimated from the data with partial sky 
coverage. The calculation of the theoretical Ci in th e con¬ 
text of a cosmological model is presented in §2.2[ This 
section also describes our modeling of the biasing rela¬ 
tion and the redshift distribution of radio galaxies. In 
0we give the maximum likelihood methodology for pa¬ 
rameter estimation. Constraints on the biasing and the 
evolution of radio activity are derived in Q In our 
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model for the redshift distribution is shown to be consis¬ 
tent with observations of the radio luminosity function 
and the source counts in the NVSS. We end in with a 
general discussion of the results. 


2. THE ANGULAR POWER SPECTRUM 

We define the surface number density contrast of 
sources 

A(f) = ■^ - 1 , (1) 

where J\f{r) is the projected number density (per stera- 
dian) in the direction f and is the mean of Af over 
the sky. In practice J\f{r) is obtained by dividing the sk y 
into small pixels (e.g. via HEALPix ( Gorski et al. 12005 )) 
and counting the number of sources in each pixel. We 
decompose A(f) in spherical harmonics Yjm as 


dim — 


L 


dnA{r)Yim(r) . 


( 2 ) 


survey 


We imagine now that dim is one out of an ensemble 
of many realizations of the same underlying continuous 


field. The variance, < \dir, 


>. 


approximated as (Peebles|[l980) 


from this ensemble is 


< \dlr. 


>.„ = 


U; “t“ I Jim ) 


(3) 


where the term involving Af (mean number density over 
the observed part of the sky) appears due to the Poisson 
discrete sampling (shot-noise) of the continuous field by 
a finite number of sources and 



Fig. 1.— The observed angular power spectrum estimated from 
the NVSS data for S > lOmJy is shown as the blue dots. The solid 
curve is the ACDM Ci obtained using the ML parameters given 
in Table 1. The dashed red curves correspond to Icj limits due to 
shot-noise and cosmic variance scatter. 


We remove 22 bright extended local radio galaxies. 


As shown by Blake & Wall (2002), the NVSS has 
significant systematic gradients m surface density 
for sources fainter than 10 mJy. Therefore, the 
analysis is restricted to sources brighter than 10 
mJy. We also discard extra bright sources {S > 
1000 mJy). 


Jim= f \Yim\^dn, (4) 

'J survey 

accounts for the partial sky coverage and gives unity for 
full sky data. The quantity Ci is the underlying angular 
power spectrum in the continuous limit (Af —>■ oo) for full 
sky coverage. It is independent of m by the assumption 
of cosmological isotropy i.e. the lack of a preferred di¬ 
rection. According to (1^, an estimate of the true power 
spectrum from the observations is 


^obs 


1 |d/mP 1 
21 + 1 ^ Jim N 


(5) 


The Icr error in this estimate is (Peebles 1980 cf. (|3| 
here) 


O’r- = 

Ci 


21 + 1 


obs I 


( 6 ) 


reflecting cosmic variance and shot-noise scatter. 


2.1. The measured Ci 

The NVSS contains ~1.7 billion sources with angular 
positions at 1.4 GHz integrated flux density S > 2.5mJy. 
The full width at half maximum resolution is 45 arc- 
sec and nearly all observations are at uniform sensitivity. 
The catalogue covers the sky north of declination —40° 
(J2000) which is almost 80% of the celestial sphere. We 
trim the data as follows ( Blake &: Wall|[2002 ). 


• Galactic sources at latitudes |5| < 5° are masked 
out to avoid Galactic contamination. 


This leaves us with 574466 sources which we use here to 
derive the observed angular power spectrum, ac¬ 

cording to ([5 1. The effect of multiple-component sources 
(e.g. radioMobes of the sam e galaxies) have be en re¬ 
moved following the recipe of Blake et al. (2004). The 
blue filled dots in figure B are for all values of I be¬ 
tween I > Iniin = 4 ancTl < Imax = 100 Corresponding 
to comoving physical scales of 50h“^ Mpc at z = 0.5. 
The partial sky coverage introduces mixing of power be¬ 
tween different Z-modes. We have assessed this mixing 
in random realizations of the surface density A, gener¬ 
ated from the AGDM cosmological model (see below) 
and found it to be negligible for I > 3. Hence, we impose 
a lower cut at I = Zmin = 4. Since we shall restrict the 
analysis to large scales, we also impose a high I threshold 
at I = Zmax = 100. In practice, for I > Imax = 100, the 
signal to noise ratio becomes very small such that no sig¬ 
nificant information is lost by the restriction to I < Umax- 
We have confirmed that in this figure agrees very 
well with Blake et al. (2004). 


2.2. Theoretical Ci 

We derive here theoretical expressions for angular 
power spectra of A(f) within a cosmological scenario. 
The redshift is denoted by z and the comoving coordinate 
by r. By z(r) we mean the redshift of a particle at a co¬ 
moving distance r in a uniform cosmological background. 
The underlying mass density held is (r, z) and the cor¬ 
responding density contrast is (5“(r, z) = p^(r, z) /p™ — 1 
where p“ is the mean background density. The density 
contrast of galaxies is denoted by 5. 
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The theoretical counterpart of A(f) is 


A(f) = / S{rr,z{r))p{r)dr , 


(7) 


where 6 is the density contrast in the three dimensional 
distribution of galaxies at a comoving coordinate r red- 
shift z = z{r). The quantity p{r)dr is the probability of 
observing a galaxy between r and (r + dr). 

The cosmological model provides the statistical proper¬ 
ties of the underlying mass density contrast field (5™(r, z), 
rather than the galaxy density contrast <5. Restricting 
the analysis to large scales (larger than a few 10 s of 
Mpc) we assume here a linear biasing relation 6{r,z) = 
b{z)6'^{r, z), as motivated by the measured three dimen- 
sion al correlation function of galaxies in redshift surveys 

e.g. Norberg et aL]|2001 Zehavi et al. 2011 Davis et al. 


20111 and analytic studies of galaxy formation (jjNusser 

et ai.|2U14). Further, on large scales, linear gravitational 
instability paradigm yields (5“(r,z) = 5J^{r)D{z) where 
D is the linear growth factor of linear fluctuati ons nor¬ 
maliz ed to unity at the present time, z = 0 (Peebles 
19801. Therefore, 


aim = j dnKYimir) ( 8 ) 

= J dnYirnir) W(r)So{rr)dr , 

where W — D{z)b{z)p{r) with z = z(r). Going to the 
Fourier domain 




(9) 


and substituting e**’’ = ji{kr)Yi%^{h)Yira{k), 

yields. 


dzm = ^ J drW J d^k5yi{kr)Yi*^{k) 

Therefore, 


Cl = < \dlr. 


> 


( 10 ) 

( 11 ) 


= ^ y dkk^Pik) 


drW ji{kr) 


and we have used < >= {2n)^d^{k—k')P{k) where 

dD is Dirac’s <5—function. The power spectrum P{k) of 
linear density fluctuations is well constrained by a va¬ 
riety of observations, primarily the cosmic mi crowave 
anisotropies measured by the Planck satellite (Planck 
Collaboration et al. 2015). The remaining ingredients 
in the theoretical niodeiing are the radial distribution 
p{r) and the biasing 6 (z). 

2.2.1. The radial distribution 


The modeling of p(r) rel ies on CENSORS )( Best et al. 


2003 Rigb y et al. 2011) and Hercules (Waddington 


et al.|| 2001 ) surveys ot radio galaxies with observed red- 


shitts and 1.4GHz fluxes. CENSORS contains 135 radio 
sources over 6 deg^ complete down to 7.2 mJy 1.4GHz 
flux limit. Spectroscopic redshifts are available for about 
73% of the sources. The remaining sources have been 



Fig. 2.— The histogram is the observed number of sources with 

5 > 7.2 mJy per redshift bin Az = 0.15. The curves “fit 1” 

6 “fit 2” are obtained using two parametric forms for the radio 
luminosity function as described in 


assigned redshifts through the I — z or K — z magnitude- 
redshift relations. In Hercules, there are 64 objects 1.2 
deg^ with S > 2 mJy. The two surveys together contain 
165 and 131 sources, respectively for S > 7.2mJy and 
10 mJy. 

The number of sources per unit redshift, N(z), in 
the two surveys jointly is the basis for our model of 
p{r) oc iV(z)dz/dr. The distribution N{z) i s represented 
by the histogram in figure According to [Rigby et al.| 
(| 2011 |), there is a difficulty in obtaining spectroscopic 
redshifts at z « 1.5, resulting in the observed gap in 
N(z). Less accurate redshift measurements compensate 
but the random error scatter objects to the sides. As we 
shall see at a later stage, the angular clustering probes 
large scale structure 0.7. Hence, the results are insen¬ 
sitive to whether galaxies near the gap are included or 
not in the modeling of p{r). Since the number of galaxies 
is too small for a robust measurement of N{z), we adopt 
the parametric form for N{z) oc p(r)dr/dz, 


iV" 


oc z ^ exp 


( 12 ) 


where the parameters oi, 02 and 03 are determined by 
maximizing the probability distribution for measuring 
N{z) as well as the observed angular power spectrum. 
Cl- Further, we impose a upper redshift cutoff in the 
abundance of sources at z = 3.5, consistent with decline 


in their number density as seen by Rigby et al. (2011). In 
this modeling we use the redshift distribution ol galaxies 
with S > lOmJy to match the flux limit of the trimmed 
version of the NVSS, as described in §2.1| The normal¬ 
ization of is such that it yields the total num¬ 

ber of sources (5 > lOmJy) in the joint CENSORS and 
Hercules data. For a Poisson distribution of N, this is 
equivalent to performing a maginalization over the nor¬ 
malization. 


2.2.2. The biasing relation 

We will resort to a variety of observational results to 
link radio galaxies to the masses of their host dark matter 
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halos. A description for the biasing of radio galaxies will 
then be obtained from well-established relations between 
the halo distribution and the mass density held. 


Best et al. (2005) studied the properties of the host 
galaxies ot a local {z < 0.1) sample of radio-loud AGN. 
They demonstrated: i) the presence of radio activity is 
a strong function, of the stellar mass, M*, of the 

host galaxies, saturating at ^ 30% for M* ^ IO^^Mq. ii) 
the shape of the radio luminosity distribution in galaxies 
with a similar stellar mass, M*, depends very we akly on 


and is close to the global lum inosity function (Mauch 
fc Sadler|2007 Best et al.|2014 ). We assume that similar 


trends extend to higher redshifts and write the fraction 
of radio loud AGNs with radio luminosity brighter than 
P as 

/,, = , (13) 


where 4>(P, 0) is close to the global radio luminosity. Fur¬ 
ther, we adopt 




ao+o-iz 


(14) 


where ao = 2.5 ± 0.2 ( Best et 5)||2005 ) will be included 
as a prior in the maximum likelihood analysis presented 
below. The fraction of radi o activity is pract ically nonex¬ 


istence at A/* < lO^°Af0 (Best et al. 2005[). One of the 


main goals of this paper is to constrain ai from the clus¬ 
tering of the NVSS galaxies and the inferred JV(z) from 
the CENSORS and Hercules. Since the stellar and halo 
masses are closely related, the dependence of radio ac¬ 
tivity on M* dictates the biasing of the radio galaxies 
through the well studied halo biasing. We assume that 
the radio luminosity depends on the halo mass, M, only 
through M* hence the form of $(P, z) is irrelevant for 
the biasing prescription. We further assume that there 
is only radio galaxy per halo. This is i ncorrect for low 
lumin osities as seen in galaxy clusters (Branches! et al. 
20051, but it is a reasonable assumption tor the lumi¬ 
nosities considered here. Therefore, we write the biasing 
factor as 


dMnh(M, z)6h(M, z)P{M,,z) 
/""dMnh(M, z)P{M,,z) 


b{z) = 


(15) 


where the stellar mass M* as a function of the halo 



later, because of the steep dependence of T on M*, the 
results are insensitive to the exact value of Mmin- At 
the high mass end we take M^ax = the mass 

scale of rich galaxy clusters. Halos with this large mass 
are rare, implying that our analysis in insensitive to this 
upper bound as well. 


Within the context of the ACDM cosmology, we seek 
constraints on the 5 parameters ao, ai, oi, 02 and 02 by 
maximizing the likelihood, Ptoti for observing (from 
NVSS) and N{z) (from the joint CENSORS and Her¬ 
cules catalogs). For brevity, we denote these parameters 
by (i = 1 • • • 5). We derive now the expression for 

Ptot- Let dim = dim/Jim ^^d continue to assume that 
the partial sky coverage does not introduce significant 
statistical correlation between dim with different {l,m). 
Restricting the analysis to large scales, the probability 
distribution function (PDF) of dim is gaussian 

(16) 
The 


(17) 

where Ci^ = {J2m d/m)/(‘^^ + !)■ Maximizing Pi with 
respect to Ci yields Ci = Ci^ — (c.f. eq. 

as expected. An estimate of the Icr error on this Ci is 
= (-d2lnP/dC/2)-i/2 = (Cz+Ar-i)V2/(2Z + l) (cf. 
Eq. 1^. We write the likelihood function for observing 
Ci‘’‘ and N{z) as 


f d/m 

2C'; + 1/A7j ’ 

where < d/^^ >= Ci + 1/N in accordance with 
PDF of the full set {dim} of all m = — P • • -I- Hs 

O 1 f 21 + 1 Ci^^ \ 

{Cl + l/A7)'+i/2 y 2 Cl+ 1 /m) 


P{dlm) 


{Cl + 1/A7)i/2 


exp 


Ptot{{e^h N, Ci^^) = P„„ Y[ P,^ n Pi , 


(18) 


i=i„ 


where Pag a (Gaussian) prior PDF incorporating the lo- 
cal, z » 0, obser vational constraint ao = 2.5 ± 0.2 of 


Best et al. (2005). Further, P^_ is a Poisson PDF for 
observing A) galaxies in the redshift bin i given a mean 
number of lV™°'^®*(2;i)Az. We will use the redshift dis¬ 
tribution of sources with S > 10mJy and Az = 0.15. 
The probability P/ is given by (|l7| ) with Ci = Ci com¬ 
puted using the |Eisenstein & Hu| (|1998|) parametrized 
form of the ACDiVl density power spectrum, P{k). The 
cosmolo gical parameters are taken from th e latest Planck 
results (Planck Gollaboration et al. 2015): Hubble con¬ 
stant Hq = 67.8kms“^Mpc“^, total matter density pa¬ 
rameter flm = 0.308, baryonic density parameter fli, = 
0.0486, linear clustering amplitude on 8h“^Mpc scale, 
tTg = 0.815 and a spectral index Ug = 0.9667. 

The evaluation of Ci with full direct numerical integra¬ 
tion is very slow. Fortunately, the Limber approximation 
is accurate to better than 5% even at small I, as revealed 
by comparison with the full calculation for some relevant 
choices of the parameters. Substitution of the Limber 
approximation 


drW{r)ji{kr) 




21 + 1 


in (111 yields, 


(19) 


3. LIKELIHOOD ANALYSIS 
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drr W'^{r) . (20) 

4. RESULTS 

We fix ^inin = 4 and ^max = 100, as the fiducial values. 
Table 1 lists the parameters 9i = obtained by max¬ 
imizing Ptot- The theoretical C; computed with these 
parameters, is shown in figure as the black solid curve. 
The dashed red lines are Icr deviations from Ci due to 
the combined error of shot-noise and cosmic variance, i.e., 
±{Ci + l/A^)-\/27(2rirT). The scatter of the blue points 
around Ci is entirely consistent with this error, as con¬ 
firmed using random samples of di^ generated assuming 
Cl is the true power spectrum. The redshift distribution 
A^modei obtained with ML parameters is plotted as the 
dotted black curve in figure This curve is remarkably 
close to the blue solid curve corresponding to a direct fit 
of to the observed redshift distribution alone. 

Conhdence levels on a single parameter 9i are obtained 
by marginalization over the remaining 4 parameters as 
follows. For a given 9i, we compute {9j} {j ^ i) which 
render a maximum in Ptot{9i, {9j:^i}). We denote by 
PiD{9i), the value Ptot evaluated at these {%}. The con¬ 
fidence limit is then Ax^{9i) = —2 \ii{Pid/P^^) where 
is the maximum of PiB{9i) which is also the max¬ 
imum Ptot with respect to all 5 parameters. The solid 
curves in the 4 panels in figure represent for ai, 
Oi, 02 and 03 . The dashed curve in each panel is an ap¬ 
proximation to PiD by means of a normal PDF. The cor¬ 
responding figure for cto turns out to be highly symmetric 
with the solid and dashed curves almost completely over¬ 
lapping and we opted to omit it here. The parameter, 
ai, the indicator for the evolution of radio activity ver¬ 
sus stellar mass, is slightly skewed to the left relative to 
a normal distribution (dashed curve). The remaining pa¬ 
rameters pertaining to are highly asymmetric. 

The super and sub scripts attached to the ML values in 
Table 1 represent the asymmetric Icr errors correspond¬ 
ing to Ax^ = 1 for the solid lines in the figure. 

The parameter ag and the corresponding error in the 
Table are very close to the values used in the prior Pag, 
indicating that ao is mainly constrained by this prior 
rather than N{z) and C°^^. Fitting iV™°'^®'(z) directly 
to the observed distribution alone without incorporating 
any information about clustering yields best fit values 
Oi = 1.199, 02 = 0.323, 03 = 0.81 in agreement with the 
ML estimates. 

The ID marginalized curves in figure do not 
contain any information on the covariance between the 
parameters. We could obtain confidence limits for 2 pa¬ 
rameters from the full Ptot by implementing a similar 
marginalization as done above for the one parameter 
case. Unfortunately, this will be extremely CPU time 
consuming. For a visual impression of the actual co- 
variance, we resort to the Metropolis-Hastings Markov 
Chain Monte Carlo (MCMC) algorithm to generate ran¬ 
dom samples of 9i drawn from Ptot{{9i}). Figure [4| plots 
7200 such random samples for the fiducial cut on / .^here 
is a strong correlation between 02 and 03 as well as oi 
and 02 . In agreement with figure the distribution of 
the individual parameters (except ao) as represented by 
the histograms deviates from normal especially for oi, 02 



and 03 . From the MCMC sampling we get < ao >= 2.52, 

< ai => 2.03, < oi >= 0.49, < 02 >= 1.21 and 

< 03 >= 1.03 with corresponding rms values of 0.19, 
0.70, 0.26, 0.45 and 0.45. These mean are within Icr de¬ 
viations from the corresponding ML estimates in Table 
1 . 

The full error matrix Cij =< {A9i){A9j) > where 
A9i = 9i — 9^^, can be estimated from the MCMC 
sampling but this requires generating an unrealistically 
large number of random sets. Hence, we approximate 
Cij as the inverse of the Hessian Hij = d^P/d9id9j at 
9i = 9f^^. The matrices Hij and Cij are listed in Table 1. 
The second row in Table 1 gives the Icr error, ^/Ca. Since 
Hij and Cij are symmetric matrices by construction the 
resulting error is also symmetric. The quantity sjCa 
represents the error when marginalization over all other 
parameters is performed. Note that, when the other pa- 
ram eter s are fixed at their ML values, the Icr scatter is 

Confidence levels in 2D pro^ctions of the parame¬ 
ter space are plotted in figure [s] as contours of Ax^ = 
—21n(Ptot/Cmax) ~ HijA9iA9j. For each pair of param¬ 
eters in each panel, Ax^ is computed by marginalizing 
over the remaining parameters. For the sake of brevity, 
we do not show projections involving 03 . This is justi¬ 
fied since it is strongly correlated with 02 and Oi as seen 
in figure The confidence levels in the ao — ai planel 
seem to agree with the distribution of points in the cor¬ 
responding panel of figure]^ The trend of decreasing ao 
with increasing ai is a reflection of the fact that a linear 
combination of these two parameters should be consis¬ 
tent with the same obse rve d clustering amplitude. 

The bias factor (Eq. 15) as a function of redshift is 
plotted in figure]^ The olas depends only on ao and ai 
and the black solid line is computed for the ML values 
of these two parameters, as indicated in the figure. The 
shaded area corresponds to ai taken within ±lcr as com¬ 
puted from the marginalized Ay^ in the top-left panel in 
figure where ao is fixed at the maximum of Ptot for 
a given ai. For comparison we also plot the case with 
no evolution in the radio activity versus stellar mass, i.e. 
(ao,Q!i) = (2.5,0). Halos of mass 3.3 x IO^^Mq (dot- 
dashed) have a bias factor close to the model prediction 
at z = 0.5, but they are less clustered at lower redshifts. 
The halo mass corresponding to the model prediction of 
b = 1.66 at z = 0 is - 4.5 X IO^^Mq. 

4.1. Dependence on l-cut and Mmin 

Different I modes of the power spectrum probe fluctu¬ 
ations at different depths (see e.g. Eq. 20). In figure 
we explore the relative contribution to Ci from fluctua¬ 
tions bet wee n r = 0 and r = r{z) i.e. with the upper 
limit in (20) replaced with r{z). The lower mode used 
in the ML, I — /min = 4, picks up about 50% of its total 
power already by z « 0.1, while / = 10 by z = 0.25. 
The modes between I = 20 and 100 gain half their power 
at redshifts between 0.35 and 0.7, with little difference 
between I = 60 and 100. In view of this, it is interest¬ 
ing to check the sensitivity of the ML estimate of ai to 
to different cuts imposed on 1. Figure shows ML es¬ 
timates of a\ for 5 different cuts on I, i.e. on /min and 
/max- Here the remaining 4 model parameters are main¬ 
tained at the ML estimates obtained with the fiducial cut 
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Fig. 3. — Confidence levels Ax^ = — 21n(PiD/P^^^) as a func¬ 
tion of individual fitting parameters (except ao). In each panel, 
marginalization is performed over the remaining parameters, as 
described in the text. 


(^min,^max) = (4,100). The horizontal bars cover these 
5 ranges in I (AZ = 20 except for the lowest point with 
Al = 16) and the vertical lines are Itr error bars. All es¬ 
timates are consistent with each other and with the ML 
value from the fiducial cut (horizontal red line) at less 
than the Icr level. There seems to be a hint for a de¬ 
creasing oi for I < 60, but it is statistically insignificant. 

So far we have assumed that radio activity exists only 
in halos with mass above = 4 x We try 

now to check whether the clustering analysis could yield 
interesting constraints on Mmin. We hx the 7 V™°dei pg^_ 
rameters, oi, 02 and 03 at the ML values in Table 1 and 
compute Ptot in tbe plane of ai and M^-m- Further, we 
marginalize over ao as usual, i.e. for a given point in 
the plane we use Oq which renders a minimum in Ptot- 
In figure we plot Ay^, defined as twice the negative 
of the logarithm of the ratio between Ptot evaluated at 
a particular point in the plane and the maximum value. 
The inner contour, corresponding to the 68 % confidence 
limit, encompasses a relatively large region, preventing 
useful constraints on the threshold halo mass. The bias 
factor, 6 , at z = 0.5 changes very little for points inside 
this confidence level provided that ao which renders a 
minimum in Ptot (i-C. the one obtained in the marginal¬ 
ization procedure) is used. This is hardly surprising since 
b has to be consistent with the observed angular power 
spectrum. 

5. CONSISTENCY WITH LUMINOSITY FUNCTION 

The joint catalog of CENSORS and Hercules used here 
for modeling N{z) contains a small number of galax¬ 
ies. It is therefore interesting to test whether our model 
Vmodei^^) is Consistent with a realistic luminosity func¬ 
tion defined as the number of sources per log(P) per co¬ 
moving Mpc^. We express the radio luminosity function, 
(j){P, z), as the sum of the st ar forming (SF) galaxie s and 


AGNs luminosity functions ( 

Mauch & Sadler||200'I 


and For the bb', we 

follow Smolcic et al. 

o 

o 

00 



Fig. 4.— Random samples drawn from Ptot (^max = 100) with 
the MCMC algorithm. 








Fig. 5.— Confidence levels in 2D projections of the parame¬ 
ter space. The contours x^ = 2.3, 6.17, and 9.21 correspond to 
68%, 95.4% and 99% confidence levels, respectively. In each panel, 
marginalization in performed over the remaining 3 parameters. 


(see also Rigby et al.||2011 1 to write 

^ P 

(l>SF {P^ Z) = 


(1 


'|2.5 


,0 


( 21 ) 


for z < Zniax = 2.5 and P*(z) = P(zmax) otherwise, and 


we use the parametrized form in |Mauch & Sadler| (|2007 1 
for ())gp at z « 0. In practice, 
small and the exact form of ())g 


the contribution of bF is 
makes little difference to 
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TABLE 1 

ML ESTIMATES OF THE MODEL 5 PARAMETERS AND THE CORRESPONDING ERROR ESTIMATES FOR /max = 100. ThE HeSSIAN (SECOND 
DERIVATIVES OF —InP WITH RESPECT TO THE FITTING PARAMETERS AT THE MLE VALUES) AND THE ERROR COVARIANCE (INVERSE OF THE 

Hessian) are also listed. 



«0 

Q;i 

ai 

Ci2 

0.3 


z.oz»_o ig4 

±0.197 

±0.742 

O-739+i;^^' 

±0.572 

0 ynt^+u.bsy 

U. ^U0_Q g4g 

±0.792 

1 ncy-l-UbDU 
i.UO/_Q gQg 

±0.531 


Hessian 


QO 

54.138 

13.431 

-80.471 

-149.24 

132.57 

Ql 

13.431 

6353 

-34.429 

-66.526 

60.32 

ai 

-80.471 

-34.429 

344.96 

610.78 

-534.96 

^2 

-149.24 

-66.526 

610.78 

1203.2 

-1125.2 

^3 

132.57 

60.32 

-534.96 

-1125.2 

1093.8 


Error covariance 


QQ 

0.039 

-0.073 

0.005 

-0.003 

-0.001 

ai 

-0.073 

0.550 

-0.129 

0.204 

0.125 

ai 

0.005 

-0.129 

0.327 

-0.444 

-0.290 

a2 

-0.003 

0.204 

-0.444 

6327 

6317 

CL3 

-0.001 

0.125 

-0.290 

0.417 

0.282 



Fig. 6.— The biasing factor of radio galaxies as a function of 
redshift. Solid curve is obtained with the ML estimates for oq 
and ai and the shaded area represents ±lcr deviation from the 
ML ai (marginalization over ckq is performed). For reference, the 
red dashed and blue dash-dot lines show, respectively, the b{z) for 
= 0 and for halos of mass 3.3 x 



z 


Fig. 7.— The fractional contribution to Ci resulting from integra¬ 
tion from redshift zero to z. The full Ci is obtained with z = 3.5. 
Only 5 values of I presented, as indicated in the figure. 



I 


Fig. 8.— The inferred ol\ from Ci for different cuts on /. The 
horizontal bar spans the I values used to obtain each point and the 
error bar represents the corresponding Icr marginalized error. 



Fig. 9.— Confidence levels in the plane of ai and the minimum 
halo mass (in solar masses) for hosting radio galaxies. 
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our analysis here. For the AGN luminosity distribution 
we write 




c 




( 22 ) 


where the redshift dependence is via C{z), P*{z), a{z) 
and j3{z) which are expanded in a Taylor series in j/ = 
ln(l + z). The coefficients of this expansion are found by 
satisfying the following constraints 

• the observed local z » 0 AGN luminosity function 
(e.g. Best et aL]|2014 ) 


• the NVSS source count, Ns{S), giving the number 
of sources per flux unit. 

• the redshift distribution of radio with S > 7.2mJy 
obtained from the 165 galaxies in Hercules and 
GENSORS. 

There is no physical or mathematical reasoning for hx- 
ing the order of the expansion for any of the redshift 
dependent functions in (/^agn- We adopt the expansions, 

P*{z) = Po + piy + P 2 y‘^ +P3y^ 

C = Co+Piy + P5y'^ + P&y^ 
a = ao+p^y + psy"^ 

P = Po- 


(23) 


Demanding that the form (22) reduces to the 
observed local luminosity function yields, C = 
10-5.33[iogioP]-lMpc-\ Po = 1024-95 WHz-1, Oq = 1.66 
and /3o = 0.42. The redshift dependence of the faint end 
slope, /?, is poorly probed by the source count, Ns{S) 
and N{z). Therefore, we simply take /3 = /?o. For a 
given choice for the parameters pi, the luminosity yields 
the number of sources per unit redshift per unit log(5') 


n{S, z) = AttDI z) 


(24) 


where Dc is the comoving distance to redshift z and 
P = P{S,z) is the luminosity corresponding to flux S 
computed for a AGDM cosmology and a power law spec¬ 
tral energy distribution, Predictions for the ob¬ 

served quantities N(z) and Ns{S) are easily derived from 
n{S,z). The parameters pi are then obtained by maxi¬ 
mizing the likelihood for observing these quantities. Ta¬ 
ble 2 gives the parameter pi along with the Hessian and 
the error covariance matrices. In “fit 1” in the table, 
galaxies with redshifts 1.45 < z < 2 were excluded in the 
N(z), while “fit 2” includes those galaxies. The AGN lu¬ 
minosity functions for these two fits computed at z = 1 


are s hown in fig ure 


from Best et al. 


(20 


10 | toge ther with other for ms taken 
"40 and Rigby et al. ( |2011[ ) (Fit 4 in 


their appendix L)j. i'here is a general agreement among 
all the curves for P > 10^®WHz“^. Figure shows the 
predictions from the luminosity functions derived here as 
the red solid (fit 1) and red dashed (fit 2) curves. There 
is a very nice agreement with the observations and the 
other approximations to N(z) plotted in the figure. An 
impressive agreement is also found between t he p redicted 
Ns{S) and the observations as seen in figure [Il| 



Fig. 10.— Various approximations to the luminosity functions. 
The red solid and dashed curves (fit 1 & fit 2) correspond to the 
luminosity functions obtained as explained in the text. 



Fig. 11. — Source count for the luminosity functions correspond¬ 
ing to fit fit 2 compared to the NVSS observations. 

6 . SUMMARY AND DISCUSSION 

We have found a mild evidence (ai > 0 at the 89% 
confidence level) for evolution in the radio activity ver¬ 
sus stellar mass relation. If all other parameters are fixed 
at their ML values, i.e. assuming that they are known we 
get a (symmetric) Icr error = 0.51. This is the 

best accuracy that the clustering of the NVSS can yield 
for ai if N(z) and ag are known. The constraint on ai is 
obtai ned by using obse rvationally motivated SHM rela¬ 
tion dMoster et al.|2013D. We have not explored the effect 
of systematic uncertaiiities in this relation which in the 
relevant redshift range (z ~ 0.3 — 0.7) r elies on a sam¬ 
ple of 28000 sources selected with Spitzer (Perez Gonzalez 


et al. 


et al. 


2008). Nonetheless, in specting hgurelO in Coupon 


(20151 shows that the Moster et al. (2013) relation 


which IS adopted here is in very good agreement with the 
majority of observations plotted in this figure. Nonethe¬ 
less we cannot ru le out the pos s ibility that ai > 0 may 
indicate that the Moster et al. (2013) SHM ratio needs 
to be slightly modihed at z ~ 0.5. Another viable inter- 
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TABLE 2 

Parameters of luminosity eunction 



Pi 

P2 

P3 

P4 

PS 

P6 

P7 

P8 

Fit 1 

0.15 

±0.14 

3.3 

±0.31 

-1.0 

±0.15 

0.34 

±0.11 

-2.33 

±0.26 

0.56 

±0.13 

-0.9 

±0.22 

1.36 

±0.32 

Fit 2 

-0.183 

±0.32 

4.32 

±0.84 

-1.635 

±0.46 

0.74 

±0.2 

-3.381 

±0.53 

1.146 

±0.29 

-0.445 

±0.24 

0.691 

±0.3 


10 Hessian (for fit 1 only) 


Pi 

0.987 

0.836 

0.806 

1.315 

1.168 

1.175 

-0.003 

0.001 

P2 

0.836 

0.716 

0.698 

1.085 

0.971 

0.985 

-0.012 

-0.005 

P3 

0.806 

0.698 

0.688 

1.019 

0.919 

0.940 

-0.020 

-0.011 

P4, 

1.315 

1.085 

1.019 

1.910 

1.664 

1.647 

0.037 

0.028 

PS 

1.168 

0.971 

0.919 

1.664 

1.459 

1.451 

0.024 

0.019 

P6 

1.175 

0.985 

0.940 

1.645 

1.451 

1.451 

0.015 

0.014 

P7 

-0.003 

-0.012 

-0.020 

0.037 

0.024 

0.015 

0.012 

0.008 

PS 

0.001 

-0.005 

-0.011 

0.028 

0.019 

0.014 

0.008 

0.005 


Error covariance (for fit 1 only) 


Pi 

0.018 

-0.041 

0.019 

-0.013 

0.028 

-0.013 

-0.004 

-0.002 

p2 

-0.041 

0.096 

-0.047 

0.027 

-0.065 

0.032 

0.022 

-0.016 

P3 

0.019 

-0.047 

0.024 

-0.012 

0.031 

-0.016 

-0.016 

0.018 

P4, 

-0.013 

0.027 

-0.012 

0.013 

-0.029 

0.014 

-0.002 

0.009 

PS 

0.028 

-0.065 

0.031 

-0.029 

0.068 

-0.034 

-0.005 

-0.007 

P6 

-0.013 

0.032 

-0.016 

0.014 

-0.034 

0.017 

0.006 

-0.003 

P7 

-0.004 

0.022 

-0.016 

-0.002 

-0.005 

0.006 

0.047 

-0.067 

PS 

-0.002 

-0.016 

0.018 

0.009 

-0.007 

-0.003 

-0.067 

0.106 


pretation for a\ > 0, is that halos hosting radio galaxies 
could be more massive than halos of other galaxies of the 
same stellar mass. This interpretation is in line with the 
conclusion of Mandelbaum et al. (20091 based on small 
scale clustering and galaxy galaxy lensing. 

The bias factor is less sensitive to the SHM since it 
has to yield a theoretical Ci which matches the observed 
angular clustering within the ACDM model. For the ML 
parameters (Table 1) we find an excellent match between 
the theoretical and observed C;. Further, estimates of a\ 
obtained for different cuts imposed on I are consistent 
with ML estimate obtained from the fiducial cut. The 
evidence for a dependence of ai on scale (cf. Fig. is 
st atistically insignificant. 


Mandelbaum et al. (20091 examined the angular clus- 
tering ot 5700 radio loud AGNs. They analyzed two 
point correlations at projected separations ;S 10fe~^ Mpc 


in the frame work of halo occupation modeling (e .g. Pea¬ 


cock & Smith 2000[ White 2001 Zheng et al. 2004) irnple- 


mented in the Millennium simulation. In contrast to our 
work which is restricted to large scales ( ^ 50/i“^ Mpc), 
their modeling depends on the division of radio galaxies 
as satellites and centrals. From the clustering measure¬ 
ment alone they obtain a mean mass of 3 x lO^^M©, for 


the central halos. The bias factor we infer (cf. Fig. 
matches that of 3.3 xJ^^Mq halos, in agreement with 


Mandelbaum et al. (2009). It is reassuring that the two 
independent studies, probing different scales, yield sim¬ 
ilar halo masses. However, incorporating galaxy-galaxy 
lensing in addition to c lustering they obtain halo masses 


(1.6 ± 0.4) X |Mandelbaum et al. (|2009|) argue 

uld 


that the deviation could be attributed to uncertainties in 
the halo modeling, which mainly affects their clustering 
results. In general, it seems that galaxy-galaxy lensing 
incorporated in halo modeling tends to push halo masses 
to lower values, as reve aled by inspecting figure 10 in 


Leauthaud et al. (2012). Nonetheless, our result is in- 
dependent of the halo model. The cause of the tension 
between the halo mass estimates based on lensing and 
clustering, remains obscure. 
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